Aharonov-Anandan phases in Lipkin-Meskov-Glick model 
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In the system of several interacting spins, geometric phases have been researched intensively. 
However, the studies are mainly focused on the adiabatic case (Berry phase), so it is necessary for 
us to study the non-adiabatic counterpart (Aharonov and Anandan phase). In this paper, we analyze 
both the non-degenerate and degenerate geometric phase of Lipkin-Meskov-Glick type model, which 
has many application in Bose-Einstein condensates and entanglement theory. Furthermore, in order 
to calculate degenerate geometric phases, the Floquet theorem and decomposition of operator are 
generalized. And the general formula is achieved. 
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I. INTRODUCTION 



Geometric phase relating to quantum mechanics is one of the most interesting developments in the recent 25 years, 
which has been discovered by Berry [lj in the context of adiabatic, unitary, cyclic evolution of time-dependent quantum 
system. He demonstrated that besides the usual dynamical phase, an additional phase relating to geometry of the 
state space was generated. Soon Simon Q give an geometrical interpretation of Berry's phase. Berry phase can be 
regarded as the holonomy of a line bundle L over the space of parameters M of the system, if L is endowed with a 
natural connection. Subsequently the degenerate case of Berry phase was generalized by Wilczek and Zee 

Discarding the assumption of adiabaticity, Aharonov and Anandan Q generalized Berry's result. The dynamical 
, ^ phase was identified as the integral of the expectation value of the Hamiltonian. The Aharonov and Anandan phase 
| (A-A phase) could be obtained by the difference between the total phase and the dynamical one and also be determined 
[ — ■ ■ by the natural connection on a U(l) principle fiber bundle over the space of projective Hilbert space. Soon after, 
>- Anandan [f| generalized the above one to the degenerate case. 

Depending on the Pancharatnam's earlier work Samuel and Bhandari 0] found a more general phase in the 
context of non-cyclic and non-unitary evolution of quantum mechanics. Furthermore, there are some reviewed papers 
@, Q and books [f| [l^l about the theoretical developments, experiments and applications of geometric phase. 

Recently, the study of geometric phase of a composite system of several spins has attracted a lot of attention. Sjovist 
fill analyzed the non-cyclic and non-adiabatic two-particle geometric phase for a pair of entangled spins in a time- 
independent uniform magnetic field. Tong, Kwek and Oh [l2j generalized the above case and calculated the geometric 
phase of the similar system in a rotating magnetic field. Yi, Wang and Zheng [l3[ investigated the adiabatic and cyclic 
geometric phase of two coupled spin- 1/2 system, one of which is driven by a varying magnetic field. Xing [l5T| 
studied further the adiabatic and cyclic geometric phase of two and three coupled spin-1/2 system with anisotropic 
interactions. Shi fill [t7| researched the Berry and A-A phases of two Heisenberg-coupled and Ising-coupled qubits 
in a rotating field respectively, both of which had a significance in quantum computing application and measuring 
geometric phases. Lately Sjovist, et. al. [Hj], analyzed the adiabatic geometric phase of ground state of finite-size 
Lipkin-Meskov-Glick type model (LMG) which consists of three spin-i particles. In this paper, we calculate the 
non-adiabatic and cyclic geometric phase, namely Aharonov and Anandan phase, of this system. 

The outline of the present paper is as follows. Section II. reviews the non-degenerate and degenerate Aharonov 
and Anandan phase. And a method of calculation for A-A phase is introduced in order to calculate the LMG 
model. Moreover, we generalize the above methods to the degenerate case. In Sec. III., the cyclic state of the 
LMG model is solved. Furthermore, the Aharonov and Anandan phase for non-degenerate and degenerate case are 
calculated respectively. In Sec. IV., the non-degenerate adiabatic phases, namely Berry phases, are analyzed. And 
the connection between A-A phase and Berry phase is discussed according to quantum adiabatic theorem. At last, a 
conclusion was drawn. 
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II. REVIEWS OF A-A PHASE AND ITS GENERALIZATION TO DEGENERACY CASE 

Consider a general quantum system with time-dependent Hamiltonian H(t), which is T- period, i.e., H(T) = H(Q). 
Given an initial state of the system ip(0), The evolution is determined by Schroinger equation, that is, 

ij t m))=H{t)\m)- (i) 

Choosing an initial state makes the evolution cyclic, i.e., ip(T) = e tx ip(0), where x is the overall phase. Furthermore, 
the overall phase can be split into two parts, namely the dynamical phase and the geometric phase [1, [H, [nj. The 
dynamical phase has a natural definition, which is 

s = - [ T (m\H(t)\m)dt. (2) 



Hence 



V = X - 6 (3) 



is the geometric phase, which is determined by purely geometric property of evolution. In order to uncover the mystery 
of the geometric phase, we introduce a single- valued vector 

m)) = e -^\m) (4) 

such that 8(0) = and 9(T) = \- Substituting Eq. (TTJ) and Eq. ^ into Eq. ^j. the dynamical can be rephrased into 
By substitution the above Equation into Eq. ([3]), one gets 



T 



Adt, (5) 



where A = (ip\-^\ip) The above formula is gauge invariant and also has a mathematical counterpart called holonomy 
in J7(l) principle fiber bundle whose base manifold is the projective Hilbert space. 

In the previous paragraph, the non-degenerate Aharonov and Anandan phase is elucidated, moreover, let's talk 
about the its generalization, which is call the degenerate A-A phase 0, As is known, the state vectors live in 
Hilbert space H. Now, we focus on the subspace V n (t) of H, where n(n > 1) is the notation of dimension of the 
subspace. Furthermore, V n (t) undergoes cyclic evolution, such that V n (T) = V n (0). Let {\ip na (t)),a = l,--- , /„} 
be an orthonormal basis of V n (t) with \ip a (T)) = |^ Q (0)) for every a. where we write \tp a (t)) instead of \tpna(t)) for 
short. By comparison with the non-degenerate case, it is easy to deduced that \ip a (t)} is degenerate generalization of 
of the single valued vector \ijj(t)). With the initial state |^ a (0)) = |f/> a (0)}, substituting 

n 

\Mt)) = ^u Pa (t)\fo(t)) 

13=1 

into the Schroinger Eq. |T]), we get 

U(t) = Texp{ [ i[A(t) - £(t)]dt}, (6) 

where A a p = ($ a (t)\4;\'4>p(t)) and £ a p = (tp a (t)\H(t)\tpp(t)) . As the matrix A and £ do not generally commute, the 
total U(T) can't be written as the product of the dynamical phase, 

r T 

fjdynamical = Texp ^_ i / S(t)dt), (7) 



(I 
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qeometric 



and the geometrical phase, 

U 9eometrtc = Texp(i I A{t)dt). (8) 
Jo 

However, the latter quantity can be transformed into the path-order integral, 

C/ 9eomet " c = Vexp(i j) A(t)dt). (9) 

Moreover, it can be regarded as the holonomy of U(N) principle fiber bundle with a natural connection, whose base 
space is the grassmann manifold. 

From the above elucidation, to calculate the Aharonov and Anandan phase is to find the cyclic initial state and 
the single-valued vector. Now we concentrate on a periodic Hamiltonian H(t) with period T. From the well-known 
Floquet theory, the time evolution operation U (t) can be written as 

u(t) = z{ty Mt , (io) 

where Z(t) is a unitary T-period operation with Z(Q) = 1 and M is a Hermitian operator. Moore and Stedman (9l.[l9| 
had discussed application of the above result into the non-degenerate case. Specifically speaking, the connection 
one-form can be expressed as 

A = i{n\Z\t)±Z(t)\n), (11) 

where |n) is the nth eigenvector of M, which is also the cyclic initial state, hence Z(t)\n) is the single valued vector. 
Moreover, we point out that this theorem can also be used to calculate the degenerate case. For simplicity, suppose 
that M has mth eigenvalue, which is degenerate and spanned up a f m subspace with orthonormal basis, and the 
the eigenvectors are expressed as \ma). It is very easy to verify that the {\ma),a = 1,2, •• • , f m } are cyclic initial 
states and single-valued vector are Z(t)\ma),a = 1,2, •■■ ,f m }- Hence, we can calculate the connection matrix of 
degenerate case, which is 

A af} = i{ma\Z\t)j t Z{t)\mP). (12) 

Thus, substituting Eq. (|lll) and Eq. (|12[) into Eq. ((SJ) and Eq. ([5]) respectively, both the non-degenerate and 
degenerate A-A phase can be obtained. 

However, the time evolution operator is very hard to get according to Eq. ©-So generally we can't follow the pro- 
cedure which is displayed in the previous paragraph to calculate the A-A phase. Nevertheless, don't be so discouraged 
and let's consider an important time-dependent Hamiltonian which has this form 

H(t) = e~ lM He lM , (13) 

where A and H are time independent. Substituting Eq. (fTB")) into Eq. ([T]),then multiplying e lAt at both side of the 
equation, one can get 

ie iAt ^-=He iAt U). 
at 

By use of the derivative formula e lAt : d\ip) / 'dt = d{e tAt \\p)) / dt — \ij})de' lAt / 'dt, it is not very difficult to see that the time 
evolution operation can be written as 

U(t) = er lAt e~ tB \ (14) 

where B = H — a saim. a gain, suppose that the Hamiltonian is H is T-period. And it is a sufficient and 
necessary condition that e~ tAT commutes with B. Hence, they have a complete set of simultaneous eigenvectors, i.e., 

B(p„ = B n 4> n , (15) 
e-^>„ = e- l9 "</)„, (16) 
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where we suppose that B have non-degenerate eigenvalues. This case had already discussed by Moore [HHl]. Its A- A 
phase |1, HI is 

t] n = (</)„| A\(f> n )T-e n . (17) 

From above, we can see that the calculational methods for the non-degenerate A-A phase was displayed. Moreover, 
there exists many degenerate A-A phase in actual physical system. Nevertheless, few physicists consider about this 
problem, except Mostafazadeh [22 |. He uses the methods of dynamical invariants to obtain the degenerate geometric 
phase. But now let us fallow another line which is more direct and convenient and generalize it to the degenerate 
case, we still confer to the non-degenerate case which is depicted above, so Eq. (fT5j) and Eq. P^)l can be expressed as 

B(j)na — b n (j) na ^ 
e -iAT, _ p-iKx 

where a = 1, 2, • • • /„ and the degenerate space is a /„-fold subspace. Subsequently, I want to transform Eq. (fT4|) into 
the formula of Eq. (|10[) . However, e~ lAt isn't generally T-periodic, so we introduce an operator ft which is commute 
with B. Hence, they have the simultaneous eigenvectors, so we can get 

Qfina = Un4>na- (18) 

Furthermore, let 51 still satisfy the following properties: 

Z(t) = e - iM e im / T , 

M = -B- fl/T. 

So the single- valued vector becomes e~ lAt e int / T \4> n a) ■ Thus we can get 

Ax0 = (4>na\A\(j) nl 3) - ^-6 a p. 

Substituting the above Eq. (fTO)) into Eq. © or Eq. ©, we can get the corresponding A-A phase 
can also be achieved by the similar calculation, i.e., 

£a/3 = {4>na\H\(j) n i3) . 

Hence, the dynamical phase can be obtained by substitution Eq. (|!?0|) into Eq. ([7]). 

III. THE AHARONOV AND ANANDAN PHASE OF LMG 

The previous section has introduced methods to calculate the non-degenerate and degenerate A-A phase. The 
methods will be use to calculate a concrete model called three qubits LMG [ll| in this section. Its Hamiltonian reads 

where 7 is an anisotropy parameter, h is the strength of an external magnetic field along the z direction, and 
S a = I Ylk=i a a ^ s t ne a th component of the spin operator(for simplicity, we set ft = 1 from now on) with o\£, <r* 
and crj are Pauli operators of the kthe qubit in the representation of a z [HI . Ignoring the trivial and constant term 
— 1(1 + 7) of H, the Hamiltonian becomes 

H = -\[o\ol + olol + o\ol + 7(4^ + a 2 a 3 + a i a 3^ --(al + a\ + a 3 z ). 

Moreover, let's consider about the isospectral one-parameter Hamiltonian family 

H = e-^'He^, (21) 



(19) 

. Another quantity 
(20) 
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whose explicit expression takes this form 



i[(cos 2 <\> + 7 sin 2 0)(er> 2 + cr 2 ^ + a\al) + (sin 2 4> + 7 cos 2 <p)(al<r 2 + er 2 er 



+ (1 — 7) sin cos < 



^4 



2 2 



^)]-|(^+^+^ 3 ) 



where <j) = tot is a varying parameter and w is the angular velocity. By a glance at Eq. (|13p . a conclusion can be 
drawn that the two Hamiltonian have a similar form. Hence, we can take the similar procedure to calculate the 
corresponding A-A phase. The time evolution operator becomes 



U(t) 



-iAt -iBt 

e e , 



(22) 



where A = ujS z and B = H — uS z . Next, we intend to calculate the cyclic initial state which is the eigen- 
vectors of B in this case. Thus, at first, the operator B must be represented in a concrete basis, which are 
{|000), |011), 1 101) , 1 110) , I 111) , 1 100), |010), |001)} , where |0) represents spin up and |1) represents spin down. So 
B takes the block-diagonal form 

u( h \ ( p (li h ^) 



where 



P(7,M) 



f-Uh + oo) -£(l- 7 ) -|(l-7) -$(l-7)\ 

4(1 -7) ^(ft + ^) -1(1 + 7) -1(1 + 7) 

4(1-7) 4(1 + 7) \{h + uj) -i(l+ 7 ) 

V -1(1-7) -1(1 + 7) -£(1+7) ^ + ^) y 



and is the 4x4 null matrix. Because B(^y,h,u>) = P(7, ft, u>) © P(7, — ft, — ui), we can investigate P(j,h,uj) and 
P(7, —ft, —cj) respectively. However, the similar information from P(7, —ft, —lu) can be obtain from P(j, ft, uf). So we 
can focus on the subspace of solutions which the sub-matrix P(7, ft, ui) acts to simplify the problem. Thus the two 
eigenvalues of P(7, ft, u) are 

Pi =-\{u + h)~{l+"t)- -Vr, 



P2 = -^(w + ft) - i(l + 7) + 3^ 

and 

P3 =P4 = \{ h + w ) + ^(1 + 7), 

where 

r = 9ft 2 + 9w 2 + 7 2 + 18ftw - 3ft7 - 37^ - 3cj - 3ft - 7 + 1. 

From above calculation, we can draw a conclusion that p\ and P2 correspond to the non-degenerate case whereas 
P3 = Pi corresponds to the degenerate case. However, its A-A phase factor is proved to be trivial. So we only focus 
on the non-degenerate case. And the corresponding two orthogonal eigenvectors are 

\<h) = -4= ( 7 + 1 - 6(w + ft) - 7 - 1 7-1 7 - 1 ) T (23) 

and 

\<h) = -4= ( 7 + 1 - 6(w + ft) + 2^ 7 - 1 7-1 7 - 1 ) T , (24) 

where 



m = 3(7 - l) 2 + [-6(w + ft) + 7 + 1 - 2v+l 2 , 



(i 



n 2 



3(7 - I) 2 + [~6(w + h) + 7 + 1 + 2<fF\ 2 , 



and T denotes the transpose operation of matrix. 

We have already calculated the cyclic initial state of the system. Moreover, from Eq. (f2"Tj) . it is very easy to verify 
that the system has 27r/o;-periodic Hamiltonian. In the following paragraphs, we will calculated the corresponding 
A-A phase of the non-degenerate case and degenerate case. 

Substituting T = 2-k/uj into e ~ tulTS ' ; by use of Eq. ([To]) , we can choose that 



Subsequently, we represent A in the given basis, which is expressed as 



.4 



/ 3 











-1 











-1 


Vo 








m 


into 


Eq 



(25) 



(26) 



the first cyclic initial state, which reads 

m = — {[1 + 7 - 6(w + h) - 

By the similar procedure, the second geometric phase is 

3tt 



(l- 7 ) 2 }-7T. 



m 



n 2 



{[1 + 7 - 6(w + h) + 2yf?f - (1 - 7 ) 2 } - 7T. 



(27) 



(28) 



The above calculation is all about A-A phase of the non-degenerate case. Next, A-A phase of the degenerate case 
will be discussed. Nevertheless, the degenerate one is trivial one, so it is necessary to reconsider another system whose 
Hamiltonian is 



H = 



According to the operator decomposition, the time evolution operator is 

U{t) = e~ lAt e- lBt , 

where A = u>S x and B = H — ujS x . Moreover, B can be represented by the given basis as 



B 



1(7-1) i( 7 -l) 1(7-1) 



|(7-1) " 4(7+1) -4(7+1) 

(7^1) -|(7 + 1) . I -1(7 + 1) 
es( ' 







[(7-1) -|(7+1) 4(7+1) 








V 















1(7-1) 

1(7-1) -|( 7 +1) 
1(7-1) -|(7+1) 



_H 

1(7-1) J(7-l) 

-g(7+l) -|(7+1) 

-I -1(7+1) 

4(7+1) -| 



Now we only focus on its two degenerate eigenvalues which are 

Bi = 1(1 + 7 - 3y/h 2 ~ 

and 



G 



B 2 = ^(1+7 + 3^ 



And the corresponding eigenvectors are 

1 



l) 



/An 



(0 ±(Vh 2 +w 2 -h) ±{h- Vh 2 + w 2 ) 0-101 



(29) 



7 



\M = -^=(o ±(VW+^-h) ±(h-V^T^) ^(VWT^-h) o -\ 1 -±) T , (30) 

V-<U2 



and 



^=(0 + ±{h+Vh 2 +lo 2 ) 0-10 l) T ; 



102 



(0 -&{h+ Vh 2 +u 2 ) ±{h + ^h 2 +oj 2 ) -±(h + Vh 2 +lo 2 ) -± 1 -i) T , 



22 



where Nu 



N 12 



N2i =4+ 4 M fe+vTO l ; and jv m = 3+ 3M^±v|!±Z). 



i+/i/V'i 2 +^ 2 

It is very easy to verify that e - l ^ G ^ ut is not a 27r/w period unitary operator, but e - t 2 ,y ^ UJt e int / T is so, where 



i o 
o 1 



(31) 



In order to calculate the degenerate connection matrix according to Eq. (|19p . we must represent A in a matrix form, 
that is 

/ 1 1 1 \ 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 
1 1 1 

V 1 1 1 / 



.4 



(32) 



Substituting Eq. ([29]), Eq. (J30J) , Eq. (J32J) and Eq. (|3Tj) into Eq. (JT9J) , the connection matrix becomes 

-1) 



2Wh 2 +uj 2 



1 
1 



(33) 



Further, by substitution Eq. (|3"3")l into Eq. (J5|), one can get the degenerate A- A phase factor 



U 



Geometric 



exp[i7r( 



Vh 2 



-1)] 



By the similar procedure, the second degenerate A-A phase factor reads 



U 



Geometric 



exp[— iir( 



Vh 2 + w 2 



1)] 



1 
1 



1 
1 



The A-A phase factor have obtained already. Furthermore, we can still get the dynamical phase factor according to 
Eq. (E0|) and Eq. 0. 



IV. DISCUSSION AND CONCLUSION 



In previous section, we have calculated the non-degenerate and degenerate A-A phase respectively. For integrity, 
it is necessary for us to concentrate on the Adiabatic phase, namely Berry phase, though Sj ovist et. al. [l8| had 
already calculated the Berry phase of the first model. But their paper probably has a minor bug of the choice of the 
instantaneous eigenvectors, which I will explain in the following paragraphs. Moreover their article only focus on the 
energetic ground state of LMG. So it is necessary for us to calculated the Berry phase for every eigenstates. But the 
degenerate Berry phase factor is still trivial, so we only concentrate on the other non-degenerate Berry phase. The 
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Hamiltonian (|2ip we focus on have a unified form expressed by Eq. (| 13|) . According Berry Wilczek and Zee 
in order to determine the Berry phase, we must calculate the instantaneous eigenstate first. The eigen equation reads 



e- iM He iM \n)=X n \n). 
It is difficult to solve the above Equation, the equation can be transformed 

H\n)' = \ n \n)\ 



(34) 



(35) 



where \n)' = e \n). For the LMG model, the above expression can be transformed to be \n) = e- % $ a * ut \n)' . In 
order to calculate Berry phase, we must choose the instantaneous eigenvectors to satisfy \n(T)) = |n.(0)). However, 
the condition can't be satisfied when the period T = 2tt/uj. After a minor modification, the eigenvector can be chosen 
to be 



(36) 



which is still satisfy Eq. (f34|) . Next, we can calculate the non-degenerate connection one-form according to Berry (I), 
which is expressed as 



A = i{n\ — \n) 



(37) 



Substituting Eq. (|36)) into Eq. (|3"7f . one can get 



Thus the Berry phase becomes 



A = iw(n|V,|»)' - J. 



T) n = I Adt. 
h 



(38) 



According to Wilczek and Zee Q, by similar calculations, the degenerate Berry phase factor can be expressed as 



U 



Geometric 



T exp(i / A(t)dt), 
Jo 



(39) 



where A a p = i{na\-i\n0) = ^uj{na\' a z \mj3)' — ^5 a /3. At first let us represent H in the given basis, which is 



H = 



( ~\h i(7~l) 1(7-1) |(7-1) \ 

J(7~l) \h -|(7 + 1) -|(7+1) 

1(7-1) -4(7+1) \h -|(7+1) 

V 1(7-1) -|(7 + 1) 4(7 + 1) \h ) 



It is not very complicated to get the eigenvalues 



Ai 



: (l+7 + 3/i + 2^), 



A 2 = --(l+7 + 3h-2V5), 



and the corresponding eigenvectors are 



ID =-L=( 1+7-6^-2^/5 7-1 7-1 7-1 ) T . 



(40) 



12} 



(1 + 7-64 + 2^7-1 7-1 7 - 1 ) T . 



TV 



(41) 
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where q = 9h 2 + 7 2 - 3/i 7 -3/i -7 + I, N t = 3( 7 - l) 2 + (1 + 7 - 6h - 2Jq) and N 2 = 3(7 - l) 2 + (1 + 7 - 6h + 2^/q). 
By the way, we ignore the degenerate eigenvalues for simplicity as well as for the triviality of the degenerate case of 
this system. Substituting Eq. (j40|) into Eq. (|38|) . we can get 



Similarly, we can also get 



^i = ^[(l+7-6/ l -2 V g) 2 -(7-l) 2 ] 



m = |f [(1 + 7 - 6h + 2^ q f - (7 - l) 2 ] - tt. 

iV2 



(42) 



(43) 



Now, let's analyze the connection between A-A phase and Berry phase further. Above all, the condition of the 
adiabatic theorem reads 



m\£H(t)\n) 



E n — E„ 



< 1. 



(44) 



Substituting Eq. ([M]) into above Eq. (|4~4")h one can simplify the above condition to 

oj \(m\a z \n)\ < 1. 

By a simple calculation, one can know |(m|<7 z |n)| ~ 1, so the condition becomes 

w < 1. 



(45) 



Substituting the above Eq. (|4"5"]) into Eq. (f!?7|) and Eq. f2"8)) , A-A phase can be reduced to Berry phase which are Eq. 
(|4"2| and Eq. (|4"3"|) . Moreover, for the degenerate A-A phase of the second model, the readers who would like to verify 
it will get the similar conclusion. 

To sum up, we have analyzed both the non-degenerate and degenerate A-A phase of LMG model. And the A-A 
phase have been reduced to corresponding Berry phase according to the condition of quantum adiabatic theorem. 
Furthermore, in order to calculate degenerate geometric phases, the Floquet theorem and decomposition of operator 
have been generalized. The general formula has been achieved. 
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Research Funds for the Central Universities. 
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